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I.  Introduction. 

Recently  there  has  been  a  revival  of  interest  in  using  Langmuir 
probes  as  a  measuring  device  in  plasmas.  However,  the  classical  low 
density  theory  cannot  be  applied  to  high  density  plasma  flows,  which 
are  encountered  in  devices  such  as  shock  tubes  and  plasma  arcs.  Ad¬ 
ditional  interest  in  the  probe  problem  arises  from  its  relation  to  the 
hypersonic  aerodynamic  problem  in  which  a  plasma  is  generated  behind 
the  bow  shock  formed  in  front  of  a  blunt  body  flying  at  hypersonic 
speeds.  In  the  probe  problem  one  is  interested  in  the  current-voltage 
characteristics  from  which  one  hopes  to  obtain  some  information  regarding 
the  properties  of  the  plasma.  In  the  blunt  body  aerodynamic  problem 
the  main  interest  is  in  the  distribution  of  charged  particles  around 
the  body  and  any  change  in  the  heat  transfer  characteristics  which  may 
occur  due  to  the  flux  of  charged  particles  to  the  body  surface  which, 
in  most  practical  situations,  is  at  the  floating  potential. 

A  continuum  theory  of  electrostatic  probes  in  a  static  isothermal 
plasma  was  given  by  Su  and  Lam''"  for  negative  probe  potentials  above  the 
floating  potential,  and  by  Cohen  for  moderate  probe  potentials  (between 
the  plasma  and  floating  potential).  In  such  analyses,  the  sheath  was 
not  assumed,  a  priori;  rather  it  turned  out  to  be  a  consequence  of  a 
careful  asymptotic  analysis.  However,  the  structure  of  the  sheath  was 
based  on  the  collision-dominated  diffusion  equation.  The  limit  of 
validity  of  such  a  description  is  obtained  by  requiring  that  the  electrical 
energy  gained  by  a  charged  particle  during  one  free  flight  is  much  less 


.■3 


-2- 


than  its  thermal  enerpy.  It  is  relatively  easy  to  show  that  such  a 
criterion  implies: 


XD  >S  1 


for  a  very  negative  probe. 


D  2/3  X 

{ — )  >>  ^ —  for  a  moderately  negative  probe  , 


where  A  ^  is  the  electron  Debye  length  based  on  the  undisturbed  charged 

particle  density,  £  is  a  typical  mean  free  path  between  charged  and 

neutral  particles,  and  r^  is  the  probe  radius.  These  inequalities  place 

rather  strong  limitations  on  the  results  in  Refs.  1  and  2.  In  general 

*  _1* 

the  ionization  fraction  has  to  be  fairly  low  (say  10  ). 

3 

These  continuum  concepts  were  later  extended  by  Lam  to  an  incom¬ 
pressible,  isothermal  flow  of  a  weakly  ionized  gas  for  moderate  surface 
potentials.  Because  of  his  assumptions  that  the  gas  is  weakly  ionized, 
incompressible,  and  isothermal,  the  diffusion  of  the  charged  particles 
to  the  solid  surface  is  decoupled  from  the  mass  motion  of  the  neutral  pas. 
The  existence  of  an  electric  field  in  the  inviscid  region  was  first  pointed 
out  in  this  work.  The  current  collected  by  the  body  is  essentially  de¬ 
termined  by  the  electron  mobility  (diffusion  due  to  the  electric  field) 

i, 

in  the  inviscid  region.  Chung  has  tackled  the  Couette  flow  and  stapna- 
tion  flow  of  weakly  ionized  pases  numerically.  The  sheath  structure  he 
obtained  checked  qualitatively  with  that  pivcn  in  Refs.  1  and  2.  This  is 
not  surprisinp  since  in  the  Couette  flow  there  is  no  convective  motion; 

* 

This  was  pointed  out  to  the  present  writer  by  J,  Fay  and  H.  F.  Frobstein 
of  ! .  Y. 
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the  poverninp  diffusion  equation  is  identical  to  the  static  case  while 
in  the  stagnation  flow  the  sheath  was  assumed  to  be  thin  nnd  close  to 
the  solid  surface  where  convective  motion  is  entirely  neplipible.  The 
analysis  in  the  inviscid  rerion  was  neglected  even  thou/rh  he  seemed  to 
note  that  there  was  residual  electric  field  intensity  at  the  outer 
edpe  of  the  viscous  boundary  layer. 

Previous  to  the  work  discussed  above,  Talbot‘S  introduced  the  concept 
of  a  "quasi-continuum"  stagnation  probe  analysis.  The  continuum  equation 
was  used  to  describe  the  diffusion  of  mass,  momentum,  and  energy  in  the 
viscous  boundary  layer,  while  within  the  sheath,  (which  was  assumed,  a 

priori,  to  occupy  a  distance  of  one  mean  free  path  from  the  surface)  the 

* 

chnrped  particles  fall  freely  down  the  potential  hill  .  The  chanpe  in 
potential  in  the  viscous  layer  was  neplected.  It  was  demonstrated  in 
l(efs.  1  and  2  that  within  the  continuum  framework  the  potential  drop  for 
a  static  plasma  outside  the  sheath  can  be  of  the  same  order  of  magnitude 
nr,  (or  lnrper  thnn)  that  within  the  sheath.  Talbot's  assumption  of  no 
chanpe  in  potential  within  the  viscous  layer  is  therefore  open  to  question. 
Kvcn  thouph  Talbot's  analysis  is  necessarily  crude,  the  quasi-continuum 
model  is  a  more  realistic  one  for  plasmus  or  hiph  ionization  fraction.  A 
complete  analysis  of  this  problem  would,  however,  require  a  kinetic  treat- 
men  t . 

* 

The  icon  of  pultinr  a  coll ision-doninnted  quasi-neutral  solution  nnd  a 
collision-five  sheath  together  was  first  suppested  in  1  yjb  by  Davydov 
and  '.,jr::n.jV:;ka.‘-'t” .  f.ucb  an  npnroxifc  nti  on  can  at  best  pive  pross  results 
such  ns  tne  current-voltnpe  characteristics. 
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In  the  present  paper,  we  shall  adopt  a  strict  continuum  description. 

The  restriction  on  such  a  continuum  analysis  which  we  have  mentioned  pre¬ 
viously  is  the  same  for  the  present  problem.  It  will  be  shown,  however, 
that  for  bodies  at  floating  potential,  the  structure  of  the  sheath  doeB 
not  enter  the  calculation  of  the  heat  transfer  and  charged  particle  dis¬ 
tribution.  Therefore  the  limitation  mentioned  is  not  relevant  to  the 
calculation  of  these  quantities.  Our  analysis  will  first  extend  Lam's 
results^  to  a  frozen  non-isothermal  plasma.  The  electric  field  in  the 
inviscid  region  will  then  be  discussed  for  a  general  flow  field.  Since 
the  flow  characteristics  in  the  viscous  boundary  layer  are  well  known, 
the  discussion  concentrates  on  the  diffusion  of  the  charged  particles  and 
the  accompanying  electric  potential  distribution.  It  will  be  shown  that 
within  the  viscous  layer,  the  diffusion  is  ambi-polar  in  nature,  even 
though  the  electron  current  is  not  necessarily  equal  to  that  for  the  ions. 

The  potential  distribution  is  decoupled  from  the  system  in  the  sense  that 
it  is  determined  after  one  has  obtained  the  solutions  for  the  other  flow 
variables.  The  probe  potential  is  assumed  to  be  moderate,  so  that  the  sheath 
is  thin  and  static  (though  with  diffusion,  of  course).  Finally  the  stagna¬ 
tion  probe  is  discussed  in  more  detail  and  an  approximate  analytic  current- 
voltage  characteristic  is  derived  under  the  assumption  of  a  very  thin  sheath. 

II.  Formulation. 

The  governing  continuum  equations  for  the  physical  system  to  be  dis¬ 


cussed  are  as  follows: 
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3v.  3v.  „  3v  3v 

11  +  3*^  "  \  fi-v  +  T~~ 

3x.  3x.  3  lk  3x  ik  3x 

xi  e  e 


q  =  -  <  grad  T  +  V  p  w  h 
«*  L.  a  -a  a 

a=l 


(perfect  gas) 


The  boundary  conditions  for  p ,  v ,  T  are  well  known.  In  addition,  we 
shall  assume  that  the  body  is  a  perfect  absorber  of  charged  particles,  so 
that  on  the  body  surface 


—  =  0  for  a  =  i  and  e. 

P 

The  potential  on  the  body  is  given  as  with  respect  to  the  potential 
far  ahead  of  the  body. 

We  shall  assume  in  our  subsequent  analysis  that  the  net  current  drawn 
by  the  body  is  small  such  that  the  Joule  heating  "£  •  E"  is  negligible. 

The  enthalpy  diffusion  flux  in  the  heat  flux  vector  £  can  be 
simplified  by  writing 

Jp  v  h  =  (h.  -  h  )  p.  w,  +  (h  -  h  )  p  w  , 

‘•a-oa  i  g  i— i  e  g'  e  —e  * 

* 

The  ion  and  electron  temperatures  T  are  assumed  to  be  equal. 
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where  the  subscript  g  stands  for  the  neutral  gas  and  h  is  the  thermal 
enthalpy  per  unit  mass.  Since 

h.  =  h  and  h  >>  h  , 

i  g  eg* 

we  have,  Yp  wh  =hp  w  =  h'  n  w,  where  h'  is  the  thermal  enthalpy 
per  electron.  We  shall  now  assume  that  the  flow  is  frozen  and  that  elec¬ 
trons  and  ions  do  not  recombine  except  on  the  body  surface.  The  enthalpy 
h  in  (4)  and  (6d)  is  then  taken  to  be  thermal  enthalpy  only.  However,  in 
our  formula  for  the  heat  transfer  to  the  wall,  h  will  include  both  the 
thermal  enthalpy  and  the  ionization  enthalpy. 

If  we  appropriately  non-dimensionalize  Eqs.  (l)  -  (5),  the  order  of 
the  magnitude  of  each  term  may  be  expressed  as  follows: 

1)  Overall  Continuity  Equation 
1:1  =  0 


3)  Momentum  Equation 


4)  Energy  Equation 
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2  2  2 
jjT  (h  +  |-)  =  |f-  +  div  {v  *  2.'  “  IT  Krad  |-*  +  “  6rad  (h  +  f~) 


-  he  {}=~  *rad  Ce  +  pa  ^a)} 
c 

p 


!  =  li  .  (1  .  i_)  I_  li  .  _i_  j  (i_  _  i) 


1  n6kT6 


Hg  ’  '*  Pr  Re  Hg  *  Pr  Re  *  'Le  VSc  Re  HjPg 


5)  Poisson  Equation 


(-i¥=  1  :  1  • 

r 

P 


Here 


Sc  =  Schmidt  number  =  (“S- )  , 

P  a  * 

D  =  Ambipolar  diffusion  coefficient  , 

pU6ro 

Re  =  Reynolds  number  =  (  —  -  * )  , 

M 

r^  =  Typical  body  dimension  , 

H  *»  lypical  total  enthalpy  , 

u.  =  Typical  velocity  , 

6 

n  *  Typical  charged  particle  number  density 

2  kT« 

A  =  Debye  length  =  --  ■■■_■  , 

D  We2 

o 


pc  b 


Le  =  Lewis  number  = 


CM 

Pr  *  Trandtl  number  =  -2— 

K 


The  reason  for  using  the  ambipolar  diffusion  coefficient,  D  ,  will  be 

a 

made  clear  in  the  discussion  of  the  viscous  boundary  layer. 
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c 
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•l 


c  c 
a  pa 


- 1 


P  c 
a  pa 

P 


i 


n  m 

c  =-e-JL 

e  p 


All  the  quantities  with  subscript  5  will  be  identified  later  as  the  values 

at  the  edge  of  the  viscous  boundary  layer. 

It  is  reasonable  to  assume  that  the  non-dimensional  parameters  Sc, 

2 

Pr,  Le,  u  g/H^,  n^kT^/p^H^  are  all  of  order  unity  in  comparison  with  the 
two  important  parameters  Re  and  X^/r^.  In  most  cases  of  practical  interest 
the  following  inequalities  are  satisfied,  i.e.. 


<< 


(7) 


Both  of  the  above  two  parameters  are  associated  with  the  relevant 
highest  order  derivatives  in  our  system  of  equations.  It  is  therefore 
expected  that  there  will  be  two  singular  perturbations  in  the  problem, 
one  for  the  viscous  layer,  associated  with  the  Reynolds  number  Re  (Prandtl 
boundary  layer)  nnd  another  for  the  sheath,  associated  with  the  Debye 
length  parameter  X^/rp  ( Langmuir  boundary  layer).  Because  of  the  in¬ 
equalities  (7),  we  see  that  the  sheath  is  imbedded  within  the  viscous  layer. 
Qualitatively,  we  can  now  say  that  there  are  three  distinct  regions  where 
different  physical  mechanisms  operate: 
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1)  Inviscid  region:  Diffusion  of  mass,  momentum,  and  energy  are 
relatively  unimportant  compared  with  convection.  Charge  neutrality  is 
maintained. 

2)  Viscous  layer:  Convection  and  diffusion  operate  simultaneously. 
Charge  neutrality  is  also  maintained  in  this  region. 

3)  Rheath:  Here  charge  separation  can  take  place.  Convection  is 
unimportant  since  the  region  i3  thin  and  adjacent  to  a  solid  surface. 
Diffusion  and  mobility  of  the  ions  and  electrons  are  the  main  features 
of  the  sheath. 

Because  of  the  assumption  that  the  sheath  region  adjacent  to  the 
surface  is  thin,  we  shall  automatically  restrict  ourselves  to  a  moderate 
probe  potential.  If  the  potential  is  strong  enough,  the  sheath  can  be 
thick  and  the  problem  of  convection  within  the  sheath  has  to  be  properly 
taken  into  account. 


III.  Inviscid  Hegion. 

In  this  region,  the  system  of  equations  becomes  doubly  degenerate. 

First  the  Lnplacian  in  the  Poisson  equation  is  neglected  on  the  basis  of 

2 

the  smallness  of  .  This  gives 


ni  3  ne  +  o  (—) 


(0) 


which  is  the  well  known  quasi-neutral  solution.  Next  we  drop  all  dissipa¬ 
tion  terms  in  Eqs.  (l)  to  (4),  i.e.. 
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+  div  (pv)  =  0  , 


p  ft  {r5  =  0  (fc}  • 


p  5t +  If =  0  (k}  • 


d  ,  1  2.  3p  ,1  . 

p  dt  (h  +  2  V  }  "  St  =  0  (FS° 


(9) 


To  this  order  of  accuracy,  the  density  p,  mass  velocity  v,  and  fluid 
enthalpy  h,  as  well  as  the  charged  particle  density  n  =  =  ng  are  de¬ 

termined  by  this  degenerate  set  of  equations.  However,  any  information 
regarding  the  electric  force  is  lost  from  the  system.  This  lost  informa¬ 
tion  can  be  recovered  by  subtracting  the  two  species  continuity  equations 
(annihilation  of  the  dominant  terms),  i.e., 


div  [n  w  -  n.  w. ]  =  0  ,  (10a) 

g  ■  g  j  * 

or  from  the  flux  relations  (Eqs.  (6a)  and  (6b) 

div  {-  p(D  -  D  )  grad  —  +  pp  (D  +  D.  )  n  grad  $}  =  0  ,  (10b) 

With  n,  p,  and  T  determined  from  Eq.  (9),  we  can  calculate  the  electric 
potential  in  the  inviscid  region  by  means  of  Eq.  (10b).  It  is  obvious  that 
the  validity  of  (10b)  is  Independent  of  the  large  Reynolds  number  assumption. 
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It  is  valid  as  long  as  quasi-neutrality  is  maintained  .  Moreover,  within 
the  sheath,  even  though  n^  ^  ng,  since  the  mass  velocity  is  small  (of  the 
order  the  ratio  of  the  sheath  thickness  to  the  viscous  layer  thickness), 

Eq.  (10b)  is  still  approximately  valid*.  Note  after  multiplying  by  e, 
the  quantity  within  the  bracket  in  (10b)  is  the  conduction  current.  Thus 
we  conclude  that  the  conduction  current  through  a  closed  surface  in  the 
flow  field  is  zero.  Within  the  viscous  and  sheath  layers  the  flux  through 
a  surface  normal  to  the  wall  is  negligible  so  that  we  have  constancy  of  the 
conduction  current  density  throughout  the  layers,  i.e., 


e  e  3  ,  T/ 
kT  p  3y  ~  ~  J/e 


(11) 


where  J  is  an  integration  constant  which  is  identified  as  the  conduction 
current  density  collected  by  the  probe,  i.e.,  J  =  =  e(r^  -  I\). 

We  have  shown  that  the  electric  potential  in  the  inviscid  region  is 
governed  by  Eq.  (10b)  with  n,  p,  and  T  obtained  from  Eqs.  (9).  One 
boundary  condition  for  Eq.  (10b)  is  obtained  by  evaluating  Eq.  (11)  at  the 
outer  edge  of  the  viscous  boundary  layer,  i.e., 

*  2 

We  have  assumed  ( \.Jr  )  <<  1/Re. 

!>  p 

^  * 

Eq.  (10b)  is  approximately  valid  within  the  sheath  if  the  convection 

there  is  negligible.  Hince  the  convective  velocity  on  the  wall  is 

zero,  the  convection  within  the  sheath  is  of  the  order  of  the  thickness 


of  the  sheath 
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kT 

r,|  'STB^H-lf-o'^-^Hl^  • 


(12a) 


In  addition,  we  require 


4>  -*■  0  at  infinity 


(12b) 


In  special  cases  such  as  stagnation  point  flow,  flow  over  a  flat 
plate,  and  the  end  wall  problem  in  a  shock  tube,  the  quantities  n,  p, 
and  T  are  approximately  constant  in  the  inviscid  region.  In  this  case 
Eq.  (10b)  and  the  boundary  condition  Eq.  (12a)  are  greatly  simplified 
so  that  the  equation  for  the  potential  and  the  boundary  conditions 
becomes 


V%  =  0  + 

(13a) 

kT 

8*  6  1  J 

ay  1  6  e(De  +  V  n 6  e  * 

(13b) 

$  -*■  0  at  infinity 

(13c) 

The  conduction  current  density  J  collected  by  the  body  is  still  an 
unknown  constant  which  must  be  determined  by  the  boundary  conditions 
specified  on  the  body  surface.  In  other  words,  J  is  determined  only 

Eq.  (10b)  is  an  elliptic  second  order  linear  partial  differential 
equation.  The  boundary  conditions  we  have  specified  uniquely  define 
n  solution. 

f 

Eq,  (13a)  was  first  shown  to  be  npproprinte  to  the  present  problem  by  Lam 
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after  one  solves  the  sheath  properly.  However,  at  the  floating  potential, 
J  is  by  definition  zero.  It  follows  then  from  (13a)  and  (13b)  that  in 
the  inviscid  region 


$  s  0 


It  is  seen  that  the  reversal  of  the  parity  of  the  electric  field  in 
the  inviscid  region  occurs  at  the  floating  potential  in  contrast  to  the 
plasma  potential  in  the  no  flow  case. 

IV.  Viscous  Layer. 

Within  the  viscous  boundary  layer,  the  dissipative  terms  (with  the 
gradient  in  the  direction  normal  to  the  wall)  are  of  the  same  order  of 
magnitude  as  the  convective  terms.  However,  quasi-neutrality  is  still  a 
good  approximation.  In  this  layer  we  have  the  following  set  of  equations: 


3p 

3t 


+  div  (pv^)  =  0  , 


Sr £>♦&£$»  . 


3  11 

a  /  e, 


3t  <r>  -  F  (“De  . 


cm 


.  du  3p  3  .  9u, 

p  dt  -  "  at  +  a7  (u  a7)  • 


r  = 0 

3y 


-15- 


From  the  above  equations  we  obtain 


1  /  O  \  3  r 

p  —  M  -  —  l 


32 - —  (rU-felp®.  !-£>] 


dt  V  3y  ll/D  +  1/D  3y  V 1  "  3y  tp  a  V 


~  [-n.w.  ]  =  |—  [-  n  w  ] 

3y  i  iy  3y  e  eyJ  * 


(15) 


-1  .  1 


where  D  =  2/(—  +  — )  is  the  ambipolnr  diffusion  coefficient, 
a  u.  U 


l  e 


V/e  conclude  from  Eq.  (15)  that 


n.w.  =  -  PD  | —  (— )  +  A 

1  iy  a  ay  P 


(16) 


n  w  =  -  pD  ^ —  (— )  +  B 
e  ey  a  3y  'p' 


(17) 


where  A  and  3  are  two  arbitrary  constants.  It  may  be  seen  from  (11)  that 
A  -  B  =  -  J/e.  It  can  also  be  seen  from  Eqs.  (l6)  and  (17)  that  the  dif¬ 
fusion  in  the  viscous  layer  is  essentially  characterized  by  ambipolnr 
diffusion.  However,  we  must  emphasize  that  the  ion  and  electron  currents 
are  the  some  only  when  the  body  is  at  the  floating  potential.  In  view  of 
Eq.  (15),  see  that  the  Schmidt  number  introduced  earlier  should  be  based 
on  the  ambipolar  diffusion  coefficient.  Since  the  latter  is  of  the  same 


Eq.  (15)  was  first  obtained  by  Chung  for  stagnation  point  flow.  However, 
he  did  not  point  out  the  simple  relations  (16)  and  (17),  which  show  that 
the  ratio  of  the  electron  and  ion  currents  cannot  be  a  constant  throughout 


the  viscous  layer 
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order  of  magnitude  as  D^,  we  conclude  that  the  thickness  of  the  diffusion 

layers  (both  for  electrons  and  ions)  are  of  the  same  order  of  magnitude 

-1/2 

as  the  viscous  momentum  layer,  i.e.,  of  order  Re  . 

The  constants  A  and  B  can  be  determined  at  the  outer  edge  of  the 
viscous  layer.  We  find  these 


B  =  -  PD  [D.e--Dl  3L  (»)  .  £.2.  ill, 
P  e  lDe  +  D.  3y  V  kT  p  3y J | 


(18) 


In  the  special  cases  such  as  stagnation  point  flow,  flat  plate  flow. 


3  ,n, 


or  the  end  wall  of  a  shock  tube,  —  (—)  i  =0.  then 

*  3y  p  |fi 


eD 


e  dip 

n,  — 


T  15 

J  e 


kTa  4  *l4  '  •  De  *  h 


(19) 


%  —  (with  D.  <<D  ) 
^  e  i  e 


Similarly 


A 


e  D  +  D. 
e  i 


J^i 

*  "  e  D 


(20) 


These  are  essentially  the  currents  due  to  the  mobilities  of  electrons 
and  ions.  It  is  seen  that  with  Dg  >>  ,  the  value  of  A  is  negligible  com¬ 

pared  with  B,  and  thus  eB  is  approximately  the  current  density  one  draws 
from  the  plasma.  At  the  floating  potential,  we  see  from  Eqs.  (19)  and 
(20)  that  A  »  B  »  0. 
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The  governing  equations  within  the  viscous  layer  can  now  he  simpli¬ 
fied  as  follows : 


■—  +  div  (pv)  =  0 


» 


P 


du  _  3£  _3_  ,  3U' 

dt  3x  +  3y  v  3y 


*£  = 
3y 


o 


2  2 
^  u  »  3p  .3  ,  , ,  1  3  ,u  «  .  u  3  ,, 

+  2  =  3t  +  3y  1_Pr  3y  2  +Pr  3y  h 


3C 

-  (7-  -  1)  pD  h  -  Bm  h  } 

Le  a  e  3y  e  e‘ 


The  electric  potential  is  obtained  from  Eq.  ( 11 )  after  this  set  of 
equations  is  solved  for  p,  n,  v,  and  T. 

Except  for  the  addition  of  the  last  term  in  the  energy  equations, 

Eqs.  (21)  are  the  usual  boundary  layer  equations  for  a  frozen  dissociated 

8  9  1 

gas  *  .  For  steady  two-dimensional  and  axi symmetric  flow  (an  analogous 

treatment  applies  to  one-dimensional  unsteady  flow) ,  we  apply  the  usual 

8  9 

similarity  transformations  used  in  compressible  boundary  layer  theory  ’  . 
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Eqs.  (23)  were  first  used  by  Talbot"*  for  the  present  type  of  problem. 
Since  he  restricted  himself  to  weakly  ionized  gases,  the  last  two  terms 
in  the  energy  equation  which  arise  from  the  energy  flux  due  to  electron 
diffusion  did  not  appear  in  his  formulation. 
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and  the  prime  indicates  differentiation  with  respect  to  n.  Eqs.  (23) 
reduce  to  the  ordinary  differential  equations  if  all  the  terms  on  the 
right  side  of  them  are  functions  of  n  alone.  For  more  detailed  discussion 
on  the  conditions  of  similarity  see  Ref.  8. 


For  a  weakly  ionized  gas  only  collisions  between  charged  and  neutral 
particles  are  taken  into  account,  so  that 


c 


_ u_ 

2pD, 


(3*/8)ptc 


For  monatomic  atoms 
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_  5  ^nnkDa  5  ^nak  *  5  o  .  2  It  8  10 

U  2  - - 2  Xtf?)—  '  3  5  “  15 

For  a  fully  ionized  gas  the  charged  particle  collisions  must  be  taken  into 
account  when  evaluating  the  mixture  transport  coefficients,  in  this  case 
one  may,  for  example,  use  the  mixture  rules  given  in  Ref.  11. 
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We  recall  that  Eqs.  (23)  are  obtained  under  the  assumption  of  quasi¬ 
neutrality.  Such  an  assumption  will  break  down  inside  the  sheath  near  the 
wall  where  we  expect  a  rapid  variation  in  potential.  Since  we  shall  later 
demonstrate  that  the  sheath  is  thin,  it  is  clear  that  the  change  in  density, 
velocity,  and  temperature  across  the  sheath  will  be  very  small.  If  we  re¬ 
quire  Eqs.  (23)  to  satisfy  the  boundary  conditions  given  on  the  wall,  the 
relative  error  introduced  in  the  solutions  for  density,  velocity,  and  tem¬ 
perature  will  be  of  order  of  the  sheath  thickness.  The  difficulty,  however, 
arises  when  we  try  to  calculate  the  potential  distribution  by  means  of  Eq. 

(11)  based  on  a  charged  particle  density  distribution  obtained  in  the 
above  fashion.  In  the  first  place,  the  current  J  cannot  be  determined. 

Second,  and  more  important,  the  electric  field  becomes  infinite  at  the  edge 
of  the  sheath.  This  suggests  that  although  the  velocity  and  temperature 
distributions  are  decoupled  from  the  system  within  the  sheath,  the  charged 
particle  density  and  potential  distribution  have  to  be  solved  simultaneously1’ 
However,  in  the  problem  of  a  blunt  body  at  floating  potential  (J  =  0),  one 
is  interested  only  in  the  charged  particle  distribution  around  the  body  and 
the  charged  particle  fluxes  to  the  wall.  As  long  as  the  sheath  is  thin, 
even  though  it  may  be  collision  dominated  or  collision  free,  the  charged 
particle  distribution  and  the  heat  flux  can  be  determined  with  an  error  of 
order  of  sheath  thickness  by  letting  Eqs.  (23)  satisfy  the  boundary  condi¬ 
tions  given  on  the  wall.  We  give  the  heat  transfer  to  a  floating  wall  as 


follows : 
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o  =  -S-l_(h+^)  +h  (1S_!jL_pd  Ii)-PD  f  (2.)  h(0) 

*v  —  3y  2  e  —  3y  H  a  3y  K  a  3y  o 


2  '  e  3y  H  a  3y 
c 

P 


.  .  n.kT 

<g'<0)  ( 


jy-  "J 
P 


p6H6 


h<°> 

( 5/2  )kT 


(2U) 


where  is  the  ionization  energy  per  electron-ion  pair.  For  the  case 

h^  >>  (5/2)kT 


n  h(0) 

C  6  6 

p 


(24a) 


The  zero  arguments  of  g  and  z  in  Eqs.  (24)  and  (24a)  are  referred  to  the 
edge  of  the  sheath,  which  in  the  floating  case  can  be  identified  as  the 
wall. 

2 

Near  the  solid  surface,  f'vp  ,  f  ''•n,  z  ^  n,  so  that  the  second 
equation  in  (23)  is  reduced  to 


(|^  z’ )’  “  0 


Thus 


z'  »  =  constant 


(25) 


or 


dn 

dy 


s  C 


(25a) 


The  corresponding  electric  potential  is  obtained  from  Eq.  (ll)  as 


It  is  not  surprising,  although  interesting,  that  the  behavior  of 

the  density  and  potential  distributions  is  exactly  the  same  as  for  the 

1  2 

quasi -neutral  solution  in  the  static  no  flow  case  *  .  The  two  constants 
C  and  D  are  related  to  the  number  density  fluxes  in  the  following  way: 


(27) 


These  constants  as  we  have  mentioned  earlier  have  to  be  determined  by 
the  boundary  conditions  at  the  wall  through  a  careful  analysis  of  the 
sheath. 

V.  Charge  Separation  Sheath. 

Since  we  shall  consider  the  case  of  moderate  potential,  the  charge 

separation  sheath  can  be  assumed  to  be  thin  .  The  analysis  within  the 

2 

sheath  will  then  be  similar  to  that  given  by  Cohen  .  The  coordinate  y 

For  a  detailed  discussion  on  this  matter,  see  the  section  on  "probes 
at  arbitrary  potentials"  in  Ref.  1. 
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in  Eqs.  (25a)  and  (26)  with  (27)  (as  is  in  the  viscous  layer)  is  the 

1/2 

physical  length  multiplied  by  Re  .  If  we  choose  the  probe  radius  r^ 

and  the  free  stream  charged  particle  density  n  to  be  the  normalization 

6 

quantities,  Eqs.  (25a)  and  (26)  with  (27)  become 


dn 

dy 


,  r  r.  r 

1  (_£.  +  _L)  _E  = 

2  ' D  D. '  n. 

e  l  6 


*e> 
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where 


j 


i.e 


r 

_E 

n. 


(28) 


(29) 
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The  electric  field  outside  the  sheath  is  of  order  Re  in  contrast 

to  order  unity  in  the  static  case.  From  Poisson's  equation,  (5),  it  can 
be  seen  that  the  stretching  factor  for  the  sheath  is  proportional  to 
(A_/r  )2^  Re1^  instead  of  (An/r  )2^  in  the  static  case.  The  electric 

Dp  Up 

field  in  the  sheath  is  of  order 


(V?/3  Re-1/3  R.1/2  -  (V2/3  Re1/6  . 

r  r 

P  P 

The  thickness  of  the  sheath  is  of  the  order 


(30) 


(V3  te1/2/R.1/J  .  (V/3  R.-1/6 

r  r 

P  P 


(31) 
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Within  the  sheath,  one  works  with  the  electric  field  instead  of 

2 

the  potential.  In  the  numerical  solutions  given  by  Cohen  ,  the  body 
potential  is  not  posed  as  the  boundary  condition;  instead  one  chooses 
a  value  of  l/2(j  +  J.),  and  the  corresponding  value  of  1/2(J  -  J, ) 

is  found  by  satisfying  the  charged  particle  density  boundary  conditions 
on  the  wall.  The  potential  on  the  body  is  then  found  by  integrating  over 
the  solution  for  the  electric  field.  In  the  present  problem  we  can  also 
assign  a  value  of  l/2(Je  +  J^),  in  which  case  l/2( is  determined 
through  the  sheath  solution  which  is  required  to  go  over  into  Eqs.  (28) 
and  (29).  With  l/2(Je  +  )  and  l/2(Je  -  J^)  known,  the  value  of  J  can 

be  calculated  and  the  electric  field  through  the  viscous  and  sheath 
layers  can  then  be  determined.  This  gives  the  potential  difference  be¬ 
tween  the  body  and  the  edge  of  the  viscous  layer.  The  potential  at  the 
latter  point,  taking  the  potential  to  be  zero  at  infinity,  is  obtained 
by  solving  Eq.  (10b)  subject  to  the  boundary  conditions  (12). 

The  probe  potential  can  be  written  as  follows : 


s  v 

where  yQ  is  the  physical  length  variable  normal  to  the  wall,  and 

1  >>  6*  >>  6  (i  =  thickness  of  sheath) 

s  s  s 

1  >>  6+  >>  6  (6  =  thickness  of  viscous  layer) 

V  V  V 
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Here  the  subscript  s  represents  the  sheath  solution;  v  the  viscous  layer 
solution;  and  I  the  inviscid  solution. 

With  a  certain  amount  of  manipulation,  the  potential  formula  can  be 
put  into  the  following  form: 


-V  in  [(-£)2/3 
P 


] 


+ 


’ 

0 


r(££)  +  rii)  _  /££)  _  i 

1  uy '  v  'dy's  ldy;v  t  s  ldyV 


dy 


(32) 


where  (^)v  -♦  g  is  the  matching  between  the  sheath  and  the  vIbcous  layer 
as  given  in  Eq,  (29).  The  logarithmic  term  in  Eq.  (32)  is  due  to  a  first 
order  pole  singularity  of  the  electric  field  in  the  matching  region  be¬ 
tween  the  sheath  and  viscous  layer.  We  display  this  logarithmic  behavior 

explicitly,  since  it  is  the  leading  term  in  Eq.  (32)  for  a  very  thin 

* 

sheath  .  In  Eq.  (32)  it  is  the  potential  at  the  outer  edge  of  the  boundary 
as  obtained  from  the  inviscid  solution. 

In  the  next  section  we  shall  apply  this  general  procedure  to  a  stag¬ 
nation  probe. 


VI .  Stagnation  Probe . 

In  the  neighborhood  of  the  stagnation  point  of  a  body,  one  can  assume 
that  the  quantities  n,  p,  and  T  are  approximately  constant  in  the  Inviscid 
region.  The  solution  of  the  electric  field  is  then  governed  by  Eq.  (13a) 
subject  to  the  boundary  conditions  (13b)  and  (13c).  Under  the  conventional 

« 

The  argument  of  the  logarithm  in  Eq.  (32)  is  the  thickness  of  the  sheath 
(see  Eq.  (31)). 
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approximation  of  stagnation  flow,  (flow  impinging  on  an  infinite  plane) 
there  is  no  solution  satisfied  by  Eq.  (13a)  subject  to  the  boundary  conr 
ditions  unless 

1)  J  =  0,  i.e.,  for  a  floating  body,  in  which  case  the  electric 
field  in  the  inviscid  region  is  identically  zero. 

2)  The  body  is  floating  except  for  a  small  but  finite  current 
element  located  at  the  stagnation  point.  Such  an  arrangement  is  of  great 
practical  interest.  In  what  follows  we  shall  consider  this  problem. 

We  shall  make  the  current  element  mentioned  above  small  enough  such 
that  the  usual  stagnation  flow  assumptions  can  be  applied  to  such  an 
element.  However,  the  element  is  considered  to  be  much  wider  than  the 
boundary  layer  thickness,  so  that  the  electric  potential  distribution  is 
essentially  one-dimensional  within  the  boundary  layer.  The  solutions  for 
the  ordinary  flow  properties  are  well  known.  From  Eqs.  (13)  the  electric 
potential  in  the  inviscid  region  is  now  given  by  the  following  equation 
and  boundary  conditions : 

V2*>  =  0 

d£  b  kT5  1_  J 
dy  *  e(De  +  D^)  n^  e 

=  0 


for  r  $  a 


y  =  0 


r  <  a 


♦  0  nt  infinity 
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where  y  =  0  is  the  outer  edge  of  the  viscous  layer,  and  where  the  current 
element  is  taken  to  be  a  circular  disc  of  radius  a.  We  are  not  very  in¬ 
terested  in  the  detailed  potential  distribution  in  the  half-space  y  i  0. 
The  only  information  required  for  the  construction  of  the  current-voltage 
characteristic  is  the  value  of  the  electric  potential  on  the  element 


y  =  0 


r  $  a 


The  solution  for  the  potential  is  given  simply  by  the  potential  dis¬ 
tribution  resulting  from  a  charged  disc  of  surface  charge  density 
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Assuming  the  variation  of  the  potential  across  the  disc  to  be  small, 
which  is  reasonable  provided  the  viscous  layer  is  thin  and  much  smaller 
than  the  dimension  of  the  current  element,  then  we  need  only  calculate 
the  potential  at  the  origin.  Ihe  potential  distribution  along  the  y-axis 
is  given  by 


$(y, 
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♦fi(y  =  0,  r  =  0)  = 


4*  A2aJ 


D  +  D 
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and 


At  the  floating  potential  $  ,  J  =  0,  so  that  $  =0  and  $  =  0  in 

f  6 

the  inviscid  region.  For  <j>  <  J  =  -  Jg  >  0,  and  by  Eq.  (33)  we 

have  <  0.  The  potential  in  the  inviscid  region  is  then  negative. 

For  potentials  slightly  above  $  ,  we  see  that  $  >  0.  The  current  in 
the  inviscid  region  is  mainly  due  to  the  electron  mobility.  Since  the 
field  is  zero  for  the  floating  potential,  the  reversal  of  parity  of  the 
electric  field  in  the  inviscid  region  occurs  at  the  floating  notential 
instead  of  the  plasma  potential  as  in  the  case  of  the  static  plasma. 
Consequently,  unlike  the  static  case,  y  =  0  is  not  a  solution  of  the 
problem,  i.e.,  when  a  probe  is  at  the  plasma  potential  there  is  still  an 
electric  field  within  the  plasma. 

Within  the  viscous  layer,  the  flow  field  is  governed  by  Eqs.  (23) 

9 

with  the  following  assumption  : 


±-  (£)  =  _£  ,  o 

dC  V<5  d e 


(3U) 


The  last  two  terms  in  the  energy  equations  of  (23)  are  due  to  the 
electron  diffusion  flux.  Under  the  stagnation  flow  assumption,  they  have 
a  similarity  property.  Ihe  solutions  for  the  velocity  and  temperature 
are  solved  in  the  usual  fashion,  i.e,,  one  ignores  the  existence  of  the 
thin  sheath  by  applying  the  wall  conditions  for  the  velocity  and 


For  a  discussion  of  this  approximation,  see  Ref.  12. 
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temperature  to  the  first  and  third  equations  of  (23).  This  procedure 
is  not  possible  for  the  second  equation  of  (23)  because  of  the  divergence 
of  the  electric  field.  For  stagnation  flow,  Eq.  (25)  is  valid  through¬ 
out  the  viscous  layer.  However  as  was  mentioned  before,  the  constants  C 
and  D  must  be  determined  by  an  analysis  of  the  sheath.  Once  the  charged 
particle  distribution  inside  the  viscous  layer  is  solved,  the  electric 
field  is  then  calculated  by  the  use  of  Eq.  (ll).  Together  with  the 
electric  field  in  the  sheath,  the  current-voltage  characteristic  1b  given 
by  Eq.  (32).  A  complete  solution  can  be  obtained  only  through  a  detailed 
numerical  computation  (both  in  the  viscous  and  the  sheath  layer).  How¬ 
ever  as  indicated  in  Eq.  (32),  the  leading  behavior  of  the  current- 
voltage  characteristic  can  be  obtained  analytically,  i.e., 

-  *  -  -  |  (Je  -  Jt)  In  [(^)2/3  Re"16)  +  0(1)  (35) 

P  P 

-U 

To  give  an  idea  of  the  accuracy  of  (35),  we  take  Ap/r  =10  , 

Re  =  10  .  The  relative  error  by  neglecting  the  terms  of  order  unity  is 
then  about  10^. 

The  potential  ^  is  normalized  by  the  thermal  energy  at  the  wall. 
From  Eqs.  (26)  to  (29),  we  have 
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Using  this  relationship  in  Eq.  (35)  we  obtain  the  current-voltage 
characteristic  as  follows: 


The  current-voltage  characteristic  is  then  given  explicitly  as 


L.  .  0.663  ^  l/2scX/3  +  ,  /r  ln  „el/6, 


enxD 
6  e 


This  formula  differs  from  the  corresponding  one  for  the  static 
plasma  first  by  the  factor  in  the  argument  of  the  logarithm. 

Furthermore,  it  is  not  difficult  to  show  for  a  spherical  probe  in  the 
static  plasma  that 

-^(Je  +  =  1 

In  the  static  plasma  we  see  then  that  the  current  collected  by  a  probe 
at  the  plasma  potential  is  the  random  current  flux,  while  in  the  flowing 
plasma,  the  current  is  given  by  the  first  term  in  Eq.  (38).  It  should  be 
pointed  out  that  Eqs.  (35)  and  (38)  become  invalid  when  the  probe 
potential  is  much  above  the  plasma  potential  such  that  “  Je  %  O’ 

VII.  Discussion. 

Within  the  framework  of  the  continuum  assumptions,  the  flow  of  an 
ionized  gas  over  a  biased  body  is  analyzed.  The  strongest  assumption  in 
the  present  work  lies  in  the  use  of  the  continuum  fluid  equations  for 
the  analysis  of  the  6heath  structure.  For  such  a  description  to  be  mean¬ 
ingful,  one  must  satisfy  the  following  condition: 
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where  l  is  the  typical  mean  free  path  in  the  system.  However,  since  the 
problem  for  a  body  at  floating  potential  is  independent  of  the  detailed 
structure  of  the  sheath,  the  above  mentioned' restriction  does  not  apply 
to  the  results  obtained  for  the  heat  flux  to  the  surface  and  for  the 
charged  particle  distribution  around  a  floating  body. 

We  have  also  assumed  that 


02  «  k  «  1 
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Under  such  an  assumption,  which  is  what  one  generally  encounters  in  the 
laboratory,  the  whole  flow  field  is  divided  into  three  regions: 

1)  Inviscid  quasi-neutral. 

2)  Viscous  (transport  of  mass,  momentum,  and  energy)  quasi-neutral. 

3)  Charge  separation  sheath.  The  body  potential  is  taken  to  be 
moderate  (between  the  plasma  and  floating  potential)  such  that  the  sheath 
is  thin  and  convection  can  be  neglected. 

The  conduction  current  is  found  to  vanish  through  any  closed  surface. 
This  is  true  in  the  outer  two  regions  because  of  the  quasi-neutrality 
(even  for  unsteady  flow).  Hie  statement  is  true  in  the  sheath  only  if 
the  convection  and  the  time  variation  inside  the  sheath  can  be  neglected. 
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As  an  example,  the  stagnation  probe  problem  is  worked  out  in  detail. 
An  explicit  expression  for  the  current-voltage  characteristic  is  given 
(Eq.  (33))  under  the  assumption  that 

2 

In  [(-r2-)2  Re1/6]  >>  1 

.As  we  have  mentioned,  an  error  of  about  10#  is  expected  in  this  formula. 

The  potential  in  the  inviscid  region  is  identically  zero  for  the 
probe  at  the  floating  potential.  Consequently,  in  contrast  to  the  static 
plasma,  the  potential  distribution  is  not  necessarily  monatomic  through- 

3 

out  the  whole  flow  field.  This  fact  was  first  pointed  out  by  Lam  . 

Also,  at  the  plasma  potential,  the  current  collected  by  the  probe  is,  in 
general,  not  the  random  flux  current. 
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